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PACS numbers: 52.20.-j, 34.80.Lx, 52.27.Gr 

In the recent work [l[ Fletcher ei aL reported on the 
novel experimental technique, enabling to measure the 
temperature of the expanding ultracold plasmas over a 
considerable time interval (up to 60-^70 /is). It was un- 
expectedly found that the electron temperature dropped 
with time as T e (t) ~ t~ a with a = 1.2 ± 0.1 s» 1 in- 
stead of a = 2, which would be expected for the adi- 
abatic cooling of electrons in the cloud expanding lin- 
early in time. The above-cited authors supposed that 
'the difference is likely due to the significant heating ef- 
fects from 3-body recombination' (i.e. the inelastic pro- 
cesses), but they did not provide sufficient quantitative 
estimates supporting such a conclusion. The aim of the 
present comment is to mention that the experimentally 
revealed i~ ^dependence can be explained under quite 
general assumptions by the purely elastic processes in the 
ultracold plasma, as it was done a few years ago in our 
work [2j. Briefly speaking, the proof of the universality 
of the t^-behavior consists of the three main steps. 

Firstly, we start from the most general form of the 
electron distribution function: 
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where r e „ and v en are the coordinates and velocities of 
electrons, and r,-„ are the coordinates of ions in a suffi- 
ciently small volume of plasma, moving with a definite 
macroscopic velocity. (The kinetic energy of ions could 
be also written here, but it is typically much less than the 
electronic one.) Despite a very complex form of the po- 
tential energy U in the regime of strong coupling, calcula- 
tion of the thermodynamic quantities depending only on 
the velocities is quite easy, because the integrals of U in 
the numerator and denominator cancel each other. Par- 
ticularly, the average kinetic energy per one particle turns 
out to be (k) = (3/2) k-oT e . This looks formally as an ex- 
pression for the ideal gas, but it is actually applicable to 
the plasma with any strength of the Coulomb's interac- 
tion U between the particles, including the state of deep 
cooling. 

Secondly, the average kinetic energy (k) can be related 
to the average potential energy (u) by the well-known 
virial relation for the Coulomb's field: (k) = (1/2) \ (u)\, 


which is also valid at the arbitrary intensity of interpar- 
ticle interactions. Strictly speaking, the virial theorem 
takes place only for the systems experiencing a finite (i.e. 
restricted in space) motion. Nevertheless, it should be 
approximately applicable also to the small (but macro- 
scopic) volume elements of the ultracold plasma cloud 
due to attractive potential by the ions, resulting in the 
quasi-confined electron motion [3]. 

Thirdly, the average potential energy can be evidently 
expressed in terms of the average interparticle distance: 
(u) ~ e 2 / (r) ~ e 2 ^ 1 / 3 , where e is the electron charge, 
and n is the concentration of charged particles. 

Finally, combining the above-written formulas, we get 
T e ~ n 1 / 3 . In particular, if the cloud expands linearly 
in time (and, consequently, its concentration changes 
as t~ 3 ), then T e ~ t~ l . Such kind of the temporal de- 
pendence was initially derived in [2J for a uniform cloud 
whose boundary moved by a linear law, but it should 
be approximately the same for a cloud with Gaussian 
density distribution whose dispersion (i.e. the character- 
istic size) increases almost linearly at large times (see 
formula (1) and below in paper [l|). 

Therefore, i~ ^dependence revealed in the above-cited 
experiment should be primarily a manifestation of "viri- 
alization" of the electron velocities in the regime of strong 
electron-ion correlations, while the heat release by inelas- 
tic processes should be of secondary importance (it might 
be responsible, particularly, for changing the exponent 
from —1 to —1.2). This is confirmed also by our ab initio 
molecular-dynamic simulations taking into account the 
strong electron-ion correlations (i.e. not based on the 
PIC method, Vlasov approximation for electrons, etc., 
which are applicable only to small-angular scattering) [3| . 
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